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ABSTRACT 

We study the red sequence in a cluster of galaxies at z = 1.62 and follow its evolution over the 
intervening 9.5 Gyr to the present day. Using deep Y JKg imaging with the HAWK-I instrument 
on the VLT we identify a tight red sequence and construct its rest-frame i-band luminosity function 
(LF). There is a marked deficit of faint red galaxies in the cluster that causes a turnover in the LF. 
We compare the red sequence LF to that for clusters at z < 0.8 correcting the luminosities for passive 
evolution. The shape of the cluster red sequence LF does not evolve between z = 1.62 and z = 0.6 
but at z < 0.6 the faint population builds up significantly. Meanwhile, between z = 1.62 to 0.6 the 
inferred total light on the red sequence grows by a factor of ~ 2 and the bright end of the LF becomes 
more populated. We construct a simple model for red sequence evolution that grows the red sequence 
in total luminosity and matches the constant LF shape at z > 0.6. In this model the cluster accretes 
quenched blue galaxies from the field and subsequently allows them to merge. We find that 3-4 
mergers among cluster galaxies during the 4 Gyr between z = 1.62 and z = 0.6 matches the observed 
luminosity function evolution between the two redshifts. The inferred merger rate is consistent with 
other studies of this cluster. Our result supports the picture that galaxy merging during the major 
growth phase of massive clusters is an important process in shaping the red sequence population at 
all luminosities. 
Subject headings: 



1. INTRODUCTION 

Understanding the formation of passive galaxies is an 
enduring problem in astronomy. These galaxies have 
very low star formation rates (SFRs), little or no cold gas, 
and dominate the population of massive g alaxies in the 
local Universe ie.s.. IKauffmann et al.ll2003l) . Their colors 
are uniformly red and they lie in a distinct region of color- 
magnitude (or color-mass) space called the red sequence. 
Studies of local passive galaxie s indicate that their stella r 
populations are very old (e.g. iBower et al.l 119921 I1998D . 
with the most massiv e passive galax ies having th e oldest 
mean stellar ages ("cglThomas et al][2 005: Gallazz Tet al.l 
[2006L iThomas et al.ll20ia but see lTrager et all (|2008| )). 
A problem with these studies is that it is not immedi- 
ately apparent how to disentangle the ages of the stars 
from the time at which they assembled into present day 
galaxies. In a dra r tiatic exam ples of the p itfalls that 
are present, iLaued ()1988[ ) and iRines et al.l (|2007[ ) find 
observational evidence for past and ongoing mergers of 
old stellar systems in brightest cluster galaxies (BCGs). 
This behavior is ec h oed i n the semi-analytic model of 
iDe Lucia fc BlaizotI (|2007[ ) who show that BCGs have 
old stellar ages but relatively recent epochs where the 
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mass was physically assembled. More recently, however, 
observations have shown that the observed stellar mass in 
BCGs has rema ined constan t since z ^ 1 (iWhilev et al.l 
[20MlCollins"et al. 2009; Sto tt et al.l[20Toll2011h . in con- 
trast to the model predictions and the implications from 
the merger remnants seen in low redshift BCGs. This 
highlights the difficulties inherent in interpreting the evo- 
lution of massive galaxies using stricly studies of the local 
universe. 

Direct lookback studies of passive galaxies shed some 
light on their or igin and evolution. For example, 
Ivan der Wei et aLl |2005) use fundamental plane obser- 
vations of field galaxies to determine that low-mass red 
galaxies at z < 1 have younger mean stellar ages than 
their more massive counterparts, similar to what is seen 
from the local studies referred to above. Efforts have also 
been made to observe the buildup of the passive popula- 
tion in situ. A p opulation of p assive galaxies is seen as 
far back as z ;^ 2 (iLabbe et al.ll20 05: Daddi ct lilll2005bl: 
'Cass ata et"all 120081 : iKriek et al.l ,2008, : ,Brammer et al.l 
2009), indicating that at least some passive galaxies 
already were in place by that time, only 3.3 Gyr af- 
ter the Big Bang. Despite their early presence on 
the stage, however, the number densitie s of passive 
galaxies evolved drama tically at z < 2 (jLabbe et al.l 
120051: IKriek et al.]|2008[) with a factor of ~ 2 in growth 
of number and mass densities at z < 1 (IBell et al 



2004; Brow n et al.] r2007: Fa ber et al.|[2007t iTavlor et al 
20091 IBrammcr ct al. 2011). A mass dependent growth 
is nonetheless debated though, as some works have 
shown that the rate of growth in the red galaxy pop- 
ulat ion at z < 1 is s lower in more massive galax- 
ies jCimatti et al.ll2Q06l: iBrown et al.ll2007t IF aber et al.l 
|2007| ) while iBrammer et al.l ( 201 If ) do not find a strong 
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mass dependence in galaxies selected to have low SFRs. 
In support of s ome mass dependence to t he growth of 
passive galaxies. iBrammer et al.l (|2011[ ) and lBundv et al.l 
I2OO6) find that there is a stellar mass at which the pas- 
sive and star- forming populations have equal density and 
that this limit evolves to lower mass with the passing of 
cosmic time. 

The most favored explanation for this growth in the 
total number of red sequence galaxies is via the transfor- 
mation of star-forming galaxies to passive ones, following 
a quenching episode (Blanton 2006; B undy et al. 200g 
iFaber et alll2007t iBrammer et~aIll201lD . In general, this 
conclusion is consistent with the observed growth in 
number density on the red sequence and the constant 
number density of blue star-forming galaxies. Unfor- 
tunately, however, the quenching mechanism has not 
been conclusively identified and there are not enough 
bright blue galaxies at < 1.5 to account for the ob- 
served evolution in the massive red sequence population, 
if these blue galaxies were sim ply to fade onto the red 
sequence fe.g iBell et al.|[200l) . One scenario to explain 
the growth of massive red galaxies without massive star- 
forming progenitors, is to grow them via the mergers of 
low-mass galaxies. These mergers have to be dissipa- 
tionless (or "dry") in order to preserve the red colors, 
isophote shapes, and low SFRs of the m assive galax- 
ies (e.g. iBell et al][200l iFaber et al.|[200l . Such a dry 
merging scenario has the added effect that it can grow 
the sizes of passive galaxies, which may explain the large 
evolution in sizes implied by some direct lo okback studies 
(iDaddi et al.ll2005at iMcIntosh et al.|[2005l: iTruiillo et a,l.l 



120061: Ivan Dokkum et al.ll2008l: Ivan der Wei et al.ll2008[) . 

From a theoretical perspective there are multiple can- 
didates for quenching, all of which involve either an ac- 
tive removal of cold gas from a galaxy, or the preven- 
tion of gas cooling onto the galaxy. Cold gas may be 
removed violently during galaxy mergers via a combina- 
tion of feedback from supernovae, gravitatio nal shocks, 
and t hat from an active galactic nucleus (e.g. iCox et all 
12004 Springclct al. 2005a). The supply of cold gas 
coming from a hot gas halo may also be shut off by 
heating from the central A GN (e.g. ICroton et aLJ 120061 : 
iMcNamara fc Nulsenll2007|) . While these processes are 
promising, it is not yet clear how well AGN can efficiently 
couple to the gas of the galaxy. Alternatively, environ- 
mental processes can either strip cold gas as a galaxy 
falls into a hot int racluster or intragrq up medium (ram 
pressure stripping; JGunn fc GottlflQT^ or strip the hot 
gas halo, thereby depriving t he galaxy of fuel fo r future 
star formation (strangulation: [Larson et al]|1980f) . These 
mechanisms for environmental quenching have been in- 
corporated into cosmological simulations by assuming 
that galaxies have their gas supply cutoff once they en- 
ter a larger dark matter halo and become a "subhalo" or 
"satellite" galaxy. The current (and simple) implemen- 
tation of this quenching has difficulty in matching the 
clustering and abundance of red galaxies, implying that 
the model is too efficie nt at quenching star formation in 
low-mass galaxies (e.g. iCoil et al.l[2008[ ). 

Clearly, an important way to constrain the ways in 
which star formation can be quenched is to study in detail 
the buildup of the red sequence as a function of galaxy 
mass or luminosity. For example, it is necessary to under- 
stand how the red sequence assembles in regimes where 



different feedback modes may dominate, e.g. as a func- 
tion of environment. Clusters of galaxies are one extreme 
of the dark matter halo power spectrum and are promis- 
ing testbeds for understanding the mechanisms by which 
star formation can be suppressed. An additional benefit 
of studying clusters is that the trip to the red sequence in 
a massive cluster is a one way street. Cooling of gas from 
the ambient intracluster medium is likely inefficient and 
the large relative velocities in massive virialized systems 
make merging unlikely (but see Ivan Dokkum et al.lll999l : 
[^an ct al. 2008i, for examples of merging red galaxies 
in a merging cluster). As we will discuss in this paper, 
however, galaxy merging may actually have been an im- 
portant player in early growth phase of clusters. 

Measurements of the luminosity function (LF) of 
red sequence galaxi e s in clusters JD e Lucia ct al. 2003 
Tanaka ct al. 2005; Stott et al. 2007; De L ucia et ail 
2007 : Gilbank et al. 2008.- Rudnick et al. 20091 ) and the 
field ([Rudnick et al.ll2009[) at z < 1 have shown unam- 
biguously that the faint red sequence population builds 
up at later times than the bright population. Because 
light correlates well with mass for red sequence galax- 
ies the luminosity trend can also be interpreted as one 
with stellar mass (although there are potential complica- 
tions with this simple M/L scaling due to the heirarchi- 
cal gr owth of galaxies on the red sequence (jSkelton et al.l 
l2011f )). This late build-up of the faint end therefore im- 
plies that whatever quenched star formation may have 
done so in low mass galaxies at preferentially later times. 
By comparing the total luminosities of clusters at 0.4 < 
z < 0.8 with their like ly descendants at z ~ from SDSS, 
iRudnick et al.l ()2009() found that the total light on the 
red sequence in clusters must have grown by a factor of 
1-3 over this span of time, si milar to t he inferred growth 
of the field red sequence (e.g. iBrown et al. 2007). 

What is still unknown from an observational stand- 
point is how the red sequence in clusters evolved at 
z > 0.8. Here we present a study of the rest-frame i- 
band red sequence and its luminosity function down to 
faint magnitudes in a z = 1.62 cluster XMM-LSS J02182- 
05102 (also known as IRC0218). This cluster was se- 
lected as an overdensity of galaxies with Spitzer /IRAC 
colors indicative of bein g at high redshi ft irrespective of 
their rest- frame colors (iPapovichj I2008D . It was subse- 
quently confirmed using spectroscopy and there are 11 
galaxies at 1.6 2 < z^rip.r. < 1-63 within 1 physical Mpc of 
the center (Pa povich et a l. 2010; Tanaka ct al. 2010). Its 
spatial structure implies that it is not relaxed and deep 
Chandra observations have found diffuse emission indica- 
tive of a potentiall y underluminous intracluster medium 
(jPierre et al.ll20lll ). This galax y cluster has a strong red 
sequence at bright magnitudes (jPapovich et al . 2010) de- 
spite not being selected on the basis of red rest-frame op- 
tical colors. In addition it has a larg e abundance of st ar- 
forming galaxies in the cluster core (|Tran et al.ll201Cll ). 

In this paper we specifically explore the red sequence 
luminosity function and its evolution. In !j2] we discuss 
the observations, and the construction of a Xj-band se- 
lected catalog. The color magnitude diagram is presented 
in fj3l We discuss the future evolutionary path of XMM- 
LSS J02182-05102 in g] The red sequence LF is pre- 
sented in [JS] along with the evolution of the LF shape and 
total luminosity of red galaxies in clusters at z < 1.62. 
We discuss our results and the implications for red se- 
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quence growth in fjHlaiid conclude in fJT) 

Throughout we assume "concordance" A-dominated 
cosmology with flu = 0.3, Oa = 0.7, and Hq = 
70 hvo km s~^ Mpc~^ unless explicitly stated otherwise. 
All magnitudes are quoted in the AB system. 

2. DATA AND OBSERVATIONS 

2.1. Image Reduction and Calibration 

We imaged XMM-LSS J02182-05 102 in the YJKs 
bands using the H AWK-I instrument (jPirard et al.ll2004l: 
iCasali et al.l 120061) on the VLT in service mode in the 
fall of 2010 (PI: Tran; 386.A-0514(A)). Our observations 
were split up into three Observing Blocks (OBs) in Y, 
two in J, and one in Kg. The total integration times 
were 9360s, 6240s, and 3000s, in YJK, respectively. The 
image quality of our observations was excellent, with 
FWHM=0'.'52, 0'.'60, and 0'.'43 in YJK, respectively. 

In order to avoid the gaps between the chips, the clus- 
ter was centered in the middle of one of the chips and the 
field of view was rotated by approximately 45 degrees to 
place other substructures on the centers of different chips. 
We dithered our observations in pseudo-random pattern 
within a box 60" on a side. 

The data were reduced using standard techniques for 
NIR imaging a nd fo ll ow th e steps taken for HAWK-I data 
from [Lidman et al.l (|2008f ) which include dark subtrac- 
tion, the application of twilight flats to remove the pixel- 
to-pixel response, sky subtraction using the XDIMSUM 
package in IRAF ^, flux calibration, astrometric calibra- 
tion, and image combination. SExtractor, SCAMP, and 
SWarp were used for the astrometric calibration and im- 
age combination *. 

The images were photometrically calibrated using an 
8" diameter aperture on 2MASS stars within the field 
for J and Kg- The Y band image was calibrated us- 
ing a UKIRT standard observed on the same night as 
our observations. The photometric zeropoints have an 
uncertainty of 0.02 mag. 

2.2. Object Detection and Photometry 

For the purposes of measuring accurate colors, the Y 
and Kg images were convolved to the FWHM of the J 
image, which had the lowest quality. This was done by 
matching the position of the stellar locus for bright, un- 
crowded, and unsaturated stars between each of the con- 
volved Y and Kg bands and the J band. 

Objects were detected from the native seeing Kg band 
image using the SExtractor software (v2.5.0). Before de- 
tection the Kg image was convolved within SExtractor 
with a Gaussian having a FWHM=5.0 pixels, or 0'.'53, 
which corresponds to the PSF size. This is applica- 
ble for the faintest sources, which are likely unresolved. 
We experimented with different combinations of the de- 
tection threshold (DETECT_THRESH) and minimum num- 
ber of required interconnected pixels (DETECT_MINAREA) 
and decided that the parameters DETECT_THRESH=2.0 
and DETECT_MINAREA=1 detected most obvious sources 
with no clear spurious detections. For the purposes of 

"^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories which are operated by the Association for Research in 
Astronomy, Inc., under the cooperative agreement with the Na- 
tional Science Foundation 

® available at http://terapix.iap.fr 



this paper, spurious sources were identified as faint Kg 
detections with no counterpart in Y or J. We chose 
DETECT_MINAREA=1 because it can be translated straight- 
forwardly into a magnitude limit for point sources, which 
simplifies understanding the detection in terms of a total 
magnitude limit (Labbe et al. 2003) . Admittedly, using 
such a low detection threshold is somewhat ambitious 
but we wanted to make sure that we went as deep as 
possible without incurring spurious counts. We therefore 
tested that most objects near our detection limit had ac- 
ceptable measurment uncertainties (see below) and that 
including the faintest objects does not affect our conclu- 
sions in any way. We quantify our detection limit below 
using completeness simulations. 

Object photometry was performed in "dual-image" 
mode, where sources were detected on the unconvolved 
Kg image and matched aperture photometry was then 
perf ormed on th e seein g-matched images. Following, 
e.g. iLabbe et al.l (|2003D , we chose "color" apertures to 
maximize the signal-to-noise (S/N) of our colors, while 
minimizing systematic errors due to crowding. For iso- 
lated sources we choose an isophotal aperture unless the 
isophotal area is less than that of a circle with d — 0.8". 
That size corresponds to 1.4xFWHM, which is the size 
that maximizes the S/N for a Gaussian profile. In the 
presence of crowding or blending we choose ad— 1" cir- 
cular aperture so that the isophotal aperture size is not 
corrupted. 

Because of correlations in the sky measurements 
caused by sub-pixel dithering, distortion corrections, and 
undetected background sources, the measured pixel-to- 
pixel rms can underestimate the true uncertainty in the 
flux of an object. Since our objects are much fainter than 
the sky, the uncertainty in the measured flux is domi- 
nated by the uncertainty in the measured value of the 
sky. Therefore, an appropriate way to estimate the flux 
uncertainty is to measure the uncertainty in the sky for 
each flux measurement. We did this using an emp t y aper - 
ture simulations described in e.g., ILabbe et al.l ()2003D . 
Briefly, for a range of aperture radii we inserted 1000 
randomly placed non-overlapping circular apertures into 
each image, excluding all objects and image boundaries. 
The distribution of measured fluxes for each empty aper- 
ture gives the real uncertainty in the sky measurement. 
As fou nd by many authors, starting with ILabbe et al.l 
(j2003| ). the measured noise scales super-linearly with 
aperture size and is significantly larger than the expecta- 
tion from pixel-to-pixel rms assuming pure uncorrelated 
Gaussian fiuctuations. For each object/band/aperture 
triplet we computed a linearized aperture size and as- 
signed the appropriate uncertainty as derived from the 
aperture simulations. Using these simulations, we mea- 
sured a formal 5-sigma limit for a d = 1" aperture of 25.2, 
24.8, and 24.1 mag in YJKg respectively, accounting for 
a point source aperture correction (see ^2.2.ip . 

We require an object to have > 50% of the effective 
exposure time to be included in the analysis. The gaps 
between the detectors is not excluded by this cut and we 
therefore only exclude regions around the edges of the 
image. 29% of the sources in the original catalog were 
excluded by this cut, but many of the sources with the 
lowest effective exposure time are likely spurious. 

Stars were rejected from the catalog using the SEx- 
tractor Stellarity index. We flagged as stars, objects with 
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Figure 1. The results of our point source completeness simulation. Left panel: The completeness fraction, defined as the ratio of detected 
sources over input sources, as a function of input magnitude. The completeness drops precipitously at Kg > 23.4 mag. Right panel: The 
error in total magnitude (including point source aperture correction) as a function of measured magnitude. The error in total magnitude 
for faint resolved sources is less than 0.1 mag. 



CLASS_STAR>0.99 or those objects with J~Ks < 0.95. 
This cut effectively removed stars on the stellar locus in 
the Ks vs. size plane. 

2.2.1. Total Magnitudes 

We measure total magnitudes in the Ks image using 
the SExtractor AUTO magnitude with a point-source 
aperture correction. The AUTO aperture scales with 
the first moment of the object radial flux profile. While 
a floor on the size of the AUTO aperture was set, for 
small (and usually faint) objects it is nonetheless the 
case that significant fiux can be missed by the small 
AUTO aperture. We compute a minimal aperture correc- 
ti qn for point sources usi ng the methodo l ogy d escribed 
in lLabbeetaD (|2003l ) and iRudnick et al.l (|2009( ). 

We found 7 bright isolated stars in our Kg image and 
calculated each of their curves of growth (COG), which 
were normalized at r = 4", corresponding to the aper- 
ture used for our zeropoint determination. These were 
then averaged to obtain an average curve of growth for 
stars in the image. For each object the aperture cor- 
rection was computed from the COG using the circu- 
larized AUTO aperture radius. For our faintest objects 
(23 < Ks{A\JTO) < 23.5) the mean aperture correction 
is 0.15 mag with an rms of 0.04 mag. This point source 
correction can be regarded as the minimal correction to 
an object's total fiux. While this correction undoubt- 
edly misses some fiux for extended objects it must be 
applied for all objects and is likely appropriate for faint 
sources near the resolution limit. Robust modeling of the 
aperture correction for the brighter and more extended 
sources is difficult as we do not have an accurate measure 
of their profile shape. For the rest of the paper, all of our 
Ks magnitudes include this aperture correction. 

2.2.2. Catalog Completeness 



To assess our catalog completeness we performed a sim- 
ulation to determine how well we detect point sources, 
which should be good analogs for the nearly unresolved 
objects near our detection limit. We added 500 sources 
to the native seeing Ks image, in batches of 100 to main- 
tain the intrinsic crowding of the images, and detected 
the objects with the same parameters used on the science 
data. We then measured the error in the magnitude as a 
function of measured magnitude and the fraction of re- 
covered sources as a function of input magnitude. The 
results of these simulations are shown in Figure [TJ Our 
90% completeness limit is Kg = 23.4 and the amount 
of flux that we miss with our total magnitude estimate 
(including aperture correction) is less than 0.1 mag for 
objects brighter than this. The 90% completeness limit 
corresponds almost exactly to the formal 5-sigma limit 
determined from our aperture simulations. It is also the 
magnitude at which the numbers of observed galaxies 
starts to rapidly fall (see Figure [2|), just as expected 
from our completeness simulation, where the complete- 
ness falls from 90 to 30% in just 0.2 mag. This rapid 
decline in completeness is a direct result of our aper- 
ture corrections, as has been seen i n other works (e.g. 
ILabbe et al.][200a IWhite et al.l[2005l) . 

3. THE COLOR MAGNITUDE DIAGRAM 

The Ks vs. Y ~ J color magnitude diagram for all 
galaxies within 2 arcminutes (1 Mpc) of the cluster center 
are show in the left panel of Figure [5] At z = 1.62 
Y ~ J corresponds to the age sensitive U — B rest-frame 
color as the Y and J filters straddle the 4000A break. 
The Ks magnitude is very close to the rest-frame i-band. 
The plotted color-magnitude diagram (CMD) includes 
the contribution of all galaxies along the line of sight. 
In constructing the LF we will statistically subtract the 
field to measure the LF of the members but here discuss 
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Figure 2. Observed color magnitude diagrams for all of the objects witliin 2 arcminutes (1 Mpc) of tlie cluster center. The catalog is 
selected in the Xg-band and the J and V-band data are 0.7 and 1.0 mag deeper than Ks, respectively. The open circles and crosses 
mark spectroscopically confirmed members and non-members respectively. Left panel The vertical blue solid line is the 90% detection 
completeness limit and the vertical dashed blue lines correspond to the formal 5 and 3cr limits in a d = 2" circular aperture. There is no 
unique relation between our i^a-band limit and a limit in y — J color but in both panels we use open squares to show those objects with a 
Y — J color uncertainty greater than 0.2 magnitudes (see right panel). The solid magenta line is a fit to all the galaxies with Y — J > 0.6 
and the dotted lines indicate ±0.3mag around this line. The labels on the top of the figure indicate the absolute i-band magnitude derived 
from passively evolving the red sequence by 1.94 magnitudes from z = 1.62 to ^ = 0. This can be directly compared with the magnitude 
in Figure |4] and is only valid for red sequence galaxies. Right panel Here the curved dashed lines indicate the reddest colors for which we 
can achieve a 3 and 5-sigma measurement on the color in a d = 0.8" circular aperture, which is appropriate for objects near our detection 
limit. As is evident from these plots, the apparent lack of red objects at faint magnitudes is real and does not stem from either detection 
incompleteness nor from an inability to measure accurate colors. 

are -0.03 ± 0.03 and 0.95 ± 0.13.^ 

We find a striking lack of red galaxies at faint mag- 
nitudes. This is reminiscient of the lack of faint red 
galaxies in c l usters at z < 1 (iDe Lucia et al.l 120041: 
Tanaka ct al. 2005; Stott ct al. 2007: 'De Lucia et al.l 
2007; Gilbank ct al. 2008; R udnick ct al. 200J) and in 
a z = 1.46 cluster ([Hilton et al. 2009), but now seen in a 
z = 1.62 cluster. It is worth addressing possible explana- 
tions for this deficit. This effect is not due to detection 
incompleteness as we detect blue galaxies down to the 
90% completeness limit. The completeness could be a 
function of color if the faint red galaxies have a larger 
size than blue galaxies at the same total magnitude but 
we find this to be an unlikely possibility given the smaller 
size of faint red g alaxies at these redshifts compare d to 
blue galaxies fe.g. lZirm et~aI]|2007l:[Toft et al.ll2007D . 

To assess the effect of color uncertainties on the per- 
ceived deficit of faint red galaxies we plot the F — J vs. 
J CMD in the right panel of Figure [51 In this plot the 
curved lines represent the reddest color to which we can 
measure the Y ~ J color at the 3 and 5-sigma level in an 
aperture of a size appropriate for objects close to our de- 
tection limit. It is important to note that this line is not 
a completeness limit but rather the limit to which col- 
ors can accurately be determined. The apparent lack of 
objects at magnitudes slightly brighter than the 5-sigma 

^ The line fit remains unchanged if we include all galaxies with 
Y — J > 0.6 and the resultant red sequence LF remains unchanged 
if we assume a constant color with magnitude for the red sequence. 



the full CMD. Even from this, however, there are several 
items worth commenting on. 

There is a clear red sequence that is well sepa- 
ra ted from the b l ue st ar-for ming galaxies , as fo und 
in iPapovich etai\ (2010) and Tanaka et al. (2010") us- 
ing significantly shallower data. Papovich et al. ( 201(| ) 
demonstrated that the brightest red sequence galaxies 
have rest-frame U — B colors and a color scatter con- 
sistent with a formation redshift of 2.35±0.1, indicating 
that they experienced their last major episode of star 
formation 1.2 ± O.lGyr before being observed. 

It is possible that some of the galaxies on the red se- 
quence have colors dominated by dust extinction since 
some of them have 24/im emission consiste nt with ob- 
scure d star formation, or possibly an AGN ([Tran et al.l 
I2010D . Yet, > 80% of the bright red sequence galaxies 
have SED fits to their rest- frame 0.15/im < A < 3.0/im 
which do not indicate significant amounts of dust ex - 
tinction ([Tran et al.l IMot IPapovichI [20llt lLot3 [2QTT1 ). 
Also, some of the spectroscopic members are on the 
red sequence and a significant fraction of these objects 
do not show evidence of emission lines in their spectra 
([Tanaka ct al. 2010), although emission li nes should have 
been detectable at the cluster redshift ([Papovich et al.l 
[20T0I) . 

We determine the slope and zeropoint of the red se- 
quence by using a robust line-fitting algorithm on all 
galaxies with Y — J > 0.6 that have errors inY — K and 
Kg less than 0.2 mag. The slope and color at Kg — 20 
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line reflects our Ks-hand detection limit and the fact 
that our J and F-band data are 0.7 and 1.0 mag deeper 
than Ks, respectively. The open symbols in both panels 
are the same objects and correspond to those with color 
errors in excess of 0.2 mag. It is clear from this plot that 
the apparent lack of objects at the faint end of the red 
sequence is neither an issue of incompleteness nor does it 
result from large color errors, as we could have measured 
colors accurately for objects in the empty faint region of 
the red sequence if they were there. 

It is important to realize that this deficit is robust to 
considerations of cluster membership as we plot all galax- 
ies within a projected distance of 1 Mpc from the cluster. 
Finally, it is possible that galaxies with large color errors 
may be preferentially scattered off the red sequence. As 
we will discuss in ij5l when performing the membership 
determination using statistical background subtraction, 
assuming all of the galaxies with large errors lie on the 
red sequence does not remove the presence of an observed 
deficit. Therefore, the lack of observed faint red galaxies 
implies that these galaxies are not present at any red- 
shift. 

4. THE FUTURE GROWTH OF XMM-LSS J02182-05102 

To place XMM-LSS J02182-05102 into an evolution- 
ary context we must identify its likely descendants at 
lower redshift. To do this we must therefore under- 
stand its expected growth in a hierarchically evolving 
Universe. There are two mass measurements of XMM- 
LSS J02182-05102. The first is derived from the x-ray 
luminosity using the local Luminosity-Temperature and 
Temperature- Mass relations, wh ich yields M ~ 7..7 ± 
3.8 X 10"M,T, ()Pierre et al.|[20lll) ^°. The second is de- 
rived from the vel ocity dispersion of th e galaxies and 
Af - 4 X IQi^'Mo ()Papovichet al.ll2010D . As shown in 
the appendix, both mass estimates are entirely consis- 
tent with that based on the total red sequence light as 
calibrated from lower redshift measurements. Given the 
large degree of substructure, this cluster is likely unre- 
laxed and we adopt the x-ray mass for the rest of the 
paper. We predict the mass growth of XM M-LSS J02182- 
05102 using the halo growth histories of iWechsler et al] 
(|2002) , which in turn uses the distrib ution of halo con- 
centrations from lBullock et ahl (l2001f) for a halo at this 
mass and at this redshift. iPoggianti et al.l ()2006l ) showed 
that these results match the predictions from the Millen- 
nium simulation (jSpringel et al.|[2005bD . In Figure [3] we 
show this growth and compare it to the masses of SDSS 
and EDisCS clusters a s derived from their g alaxy veloc- 
ity dispersions (Milvang- Jensen et al.l[2008[ ). It is clear 
that XMM-LSS J02182-05102 will evolve into a typical 
(log(M/MQ) ~ 14.3 - 15) system by z < 1 even tak- 
ing into account the dispersion in merger histories of 
such massive halos. We will use this projected growth 
in subsequent sections when studying the evolution of 
the galaxies in XMM-LSS J02182-05102." 

^^ The x-ray luminosity for this cluster is not given in 
IPierre et all i20rih . 

^^ Had we adopted the mass from the galaxy velocity dispersion, 
the predicted descendants of XMM-LSS J02182-05102 would have 
been among the more massive clusters in the Universe at any red- 
shift. On one hand, such a massive nature for XMM-LSS J02182- 
05102 may be unlikely given how rare such objects are. Also, it 
is worth noting that the mass function of clusters at any epoch is 



tlookbaok/[Gyr] 

13 5 7 



10 



16.5 



16 



15.5 



I I I 

I XMM-LSS J02182- 
hA EDisCS clusters 
■ SDSS clusters 



05102 



O 




Figure 3. The expected growth in mass of XMM-LSS J02182- 
05102 to z = 0. The red solid square repr esents the x -ray derived 
mass for XMM-LSS J02182-05102 dPierre et aLll^OTiD . The open 
square connected to the soli d by a dotted l i ne ind icates the dy- 
namical mass estimate from IPapovich et al.l I I201CI) , which is sig- 
nificantly higher than the x-ray mass. The disagreement between 
the two likely reflects the unrelaxed dynamical state of the system. 
The gray triangles represent the EDisCS clusters. The shaded his- 
togram represents the distribution of dynamical masses from our 
local SDSS sample. The diagonal dashed line gives the median ex- 
pected growth in mass for XMM-LSS J02182-05102 from the halo 
growth histories of Wechslcr ct al. ( 2002) and using the halo con- 
centrations of Bullock ct al. (2001). The dotted lines and yellow 
shaded region give the 68% confidence limits on the growth given 
by the range in the halo concentrations for a cluster of this mass 
observed a.t z = 1.62. The likely descendants of XMM-LSS J02182- 
05102 will be fairly typical clusters at all redshifts. 

5. THE RED SEQUENCE LUMINOSITY FUNCTION 

5.L Cluster Membership 

As very few spectroscopic redshifts are available in the 
core of the cluster, and most of them are for blue galax- 
ies, we identify which galaxies are red seq uence members 
using a statistical subtraction technique. iRudnick et al.l 
([2009) demonstrated that using statistical background 
subtraction to isolate red sequence members in clusters 
at 0.4 < z < 0.8 yielded results that were identical 
to those computed using accurate photometric redshifts. 
An advantage of statistical background subtraction for 
membership of faint galaxies is that photometric redshift- 
based membership techniques often rely on integrating 
the redshift probability distribution (P(z)). Faint galax- 
ies will have lower photometric S/N and hence broader 
P{z) and the traditional method of establishing a thresh- 
old in the integrated probability (jBrunner fc Lubiiill2000f) 
will more likely reject these object, which is undesirable 
given the goals of this paper. On the other hand, sta- 
tistical background subtraction does not have a formal 

very steep and that low-mass descendants will be preferred over 
high-mass ones. Therefore, in addition to more precise mass mea- 
surements, a more careful treatment of the cluster mass function 
would be required to obtain a more accurate estimate of the typical 
descendant at each redshift. 
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Table 1 

Rest-frame i-band Luminosity Function of 
XMM-LSS J02182-05102 
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^i,bright 


Mijaiut 


# 


5<I>_ 


5<S>+ 


-24.5 


-24.0 


1.0 


0.83 


2.41 


-24.0 


-23.5 


0.85 


0.85 


2.36 


-23.5 


-23.0 


5.08 


2.17 


3.47 


-23.0 


-22.5 


5.85 


2.34 


3.62 


-22.5 


-22.0 


2.79 


1.56 


2.94 


-22.0 


-21.5 








2.0 


-21.5 


-21.0 


2.19 


1.36 


2.78 


-21.0 


-20.5 


0.54 


0.54 


2.24 



Note. — <I> corresponds to the number 
of galaxies in the magnitude range speci- 
fied with no evolution corrections applied. 
The non-integer values are a byproduct of 
our statistical background subtraction tech- 
nique. The last two columns are the positive 
and negative uncertainty in this number as 
determined from Poisson uncertainty on the 
counts. There is one magnitude bin with no 
galaxies detected. 

dependence on the S/N, other than being optimal when 
the distribution of color and magnitude errors of the field 
and cluster samples are identical, which is the case for 
our analysis. 

Based on our experience () Rudnick et al .1120091) a suit- 
able background subtraction field needs to have compa- 
rable bandpasses (for red sequence selection) , depth, and 
total magnitude measurements to the main survey field. 
A wide area is also desireable. We searched the litera- 
ture for data that satisfy these requirements but have not 
found them. For example, UKIDSS and GOODS while 
deep and wide, do not include Y-band data at our depth. 
We therefore decide to use the outskirts of our HAWK-I 
image to define the background population. While it is 
possible that the outskirts of our image may contain the 
imprint of the associated large scale structure, it may in- 
deed be correct to subtract this "local" field rather than 
a field drawn from a cosmic average. 

For the purposes of the subtraction, we define two re- 
gions, the cluster at rproj < 0.75 Mpc^^ and the field at 
fproj > 1.5 Mpc. After selecting galaxies within ±0.3mag 
of the red sequence, we bin the observed Kg band mag- 
nitudes in both regions and then subtract the field his- 
togram from the cluster histogram, normalizing the field 
histogram by the ratio of the field area to the cluster 
area. This results in an observed Kg LF for likely red 
sequence cluster members. 

5.2. The Luminosity Function 

We construct the LF in the rest-frame z -band for com- 
parison with lower redshift studies (e.g. IRudnick et al.l 
[2009) . To do this we take advantage of the fact that the 
Kg band is very close to the redshifted rest-frame z -band 
filter. Using the technique described in Rudnic k et al.l 
([2003) to compute irest we find that the k-corrections 
for red sequence galaxies at the cluster redshift are con- 
sistent with a constant tuk^ — Mi = 44.13 with a 0.03 
mag dispersion. We apply this correction to the observed 
LF to obtain a rest-frame Vest-band LF for red sequence 

^^ This is chosen to match the aperture used for the lower red- 
shift comparison samples. Its exact value does not affect our con- 
clusions. 
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Figure 4. The rest-frame i-band luminosity function for the clus- 
ter XMM-LSS J02182-05102 at z < 1.62 is shown as the solid red 
circles. Magnitude bins with no objects are not plotted. The gray 
open circles, squares, and triangle symbols are the LFs for com- 
posites of cluster s at z < 0.06, . 4 < z < 0.6, and 0.6 < z < 0.8 
respectively from lKiidnick et al.l II2009I ). The error bars are Pois- 
son errors only. The dotted, dash, and dot-dashed lines are the 
ISchechteil II1976I ) fits to the z < 0.8 clusters with a left as a free 
parameter. The solid line is the fit to the z = 1.62 cluster with a 
left fixed to the z = 0.7 value. The magnitudes shown here have all 
been passively evolved to z = as described in the text. The LFs 
have also been scaled vertically to have the same integrated lumi- 
nosity. The measured LF at 2 = 1.62 is similar in shape to that at 
z = 0.7 and shows a significant decline towards fainter magnitudes 
and a lack of bright galaxies. 

cluster members. 

The shape of this LF is insensitive to the exact slope of 
the red sequence and the width of our red sequence cut. 
We also test how the LF would change if we assumed that 
all galaxies with Y — J errors more than 0.2 mag were all 
on the red sequence. Since this adds similar numbers of 
galaxies to the field and cluster CMDs the shape of the 
subtracted cluster LF is unchanged. 

In Table [1] we give the red sequence LF of XMM- 
LSS J02182-05102 and show it in Figure [1 compared to 
the composite red sequence luminosity fu nction in clus- 
ters a t z = 0, z = 0.5, and z = 0.7 from IRudnick et all 
(|2009f ). As we will discuss shortly, there is clearly strong 
evolution in the shape of the LF at z < 0.7, with little 
shape evolution from z = 1.62 to z = 0.7. 

In iSJwc showed that XMM-LSS J02182-05102 is hkely 
the progenitor of "typical" clusters in EDisCS and SDSS 
and so in Figure [H we therefore compare the XMM- 
LSS J02182-05102 LF to the total composites for these 
two lower redshift samples. In Figure [4| we have faded 
all the LFs to z = assuming that the red sequence 
galaxies all formed at z = 2.35 and evolved passively 
thereafter. As has been seen before, the red sequence 
LF evolves strongly in shape at z < 0.8 with the 



faint galaxy po pulation building up progressively to wards 
lower r edshift ('De Lucia et al.ll20q4 iT anaka et al."2005l: 
iStott et al., ,2007; De Luci a et al.l 12007; .Gilbank et aJl 
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I2008t iRudnicket alll2009[ ). This has been associated 
with the late addition of faint galaxies to the red se- 
quence. 

The assumed fading is consistent with the evolution 
in the mean color and color s catter of the red se - 
quence in this a n d oth er clusters (jPapovich et al.ir201Cll ). 
iWhitaker et al.l (|2Q1CI[ ) found that passive galaxies at 
z ~ 1.6 have a spread in their ages of about 1 Gyr but 
this would result in a spread in the measured fading to 
z = of 0.6 mag and should not affect our results. The 
spread in fading is even less if the galaxies are only faded 
to z = 0.7 as we do in §5.31 

Given this strong evolution at z < 0.8 it is striking 
that the observed LF at z = 1.62 is remarkably similar 
in shape to that z = 0.7, with the same turnover to- 
wards faint magnitudes. At the same time, the bright 
end of the LF at z = 1.62 appears to be underpopu- 
lated with respect to the individual lower redshift cluster s 
(see individual cluster LFs from iRudnick et all I20090 . 
This is similar to what h as been found in a z = 1.46 
cluster [Hilton et all ()20Q9[ ). Taken together this implies 
that galaxies must be added to the bright end between 
z = 1.62 and z = 0.7 but that the shape of the LF should 
remain relatively constant. In the following section we 
will discuss possible scenarios for the evolution of the 
LF and how it relates to the evolution of the cluster red 
sequence. 

The iSchechted (|1976l ) funct ion fits to the EDisCS and 
SDSS LFs were presented in Rudnick ct al. (2009) and 
were computed with a, M* and the normalization as free 
parameters. Our LF at z = 1.62 does not have enough 
signal-to-noise to allow a simultaneous determination of 
a and M* but we attempted to constrain M* by fixing 
a = 0.17 as determined from the z = 0.7 clusters. While 
the lack of bright galaxies allowed no strong constraints 
to be placed on M*, the best-fit Schechter function with 
a fixed a = 0.17 (Figure [3]) are statistically acceptable 
fi ts to the observe LF 

iTanaka et al.l (|20100 showed that the red sequence LF 
of XMM-LSS J02182-05102 as derived from significantly 
shallower NIR data appeared simila r to that from groups 
at z = 1.1 from lTanaka et aP ()2008[ ) which in turn had a 
deficit of faint galaxies. At the same time, the very mas- 
sive z = 1.1 cluster from iTanaka et all (|2008f ) had a red 
sequence LF that appeared similar to that from SDSS 
clusters, but was very different from that for XMM- 
LSS J02182-05102. The strong cluster m ass dependence 
of the LF contrasts to the result from IRudnick et al.l 
(|2009f ) who found only a very weak dependence of red 
sequence LF shape on cluster mass. Indeed, the LF 
for the most massive cluster from EDisCS in its high 
redshift sample (CL1216. 8-1201, z — 0.8) is consistent 
with both the z — 0.7 EDisCS composite LF and that 
for XMM-LSS J02182-05102 but inconsis tent with the 
z = .5 E DisCS composite. As noted in ITanaka et al.l 
(|2008D and IRudnick et al.l ()2009f ). this apparent discrep- 
ancy may result from the large cluster-to-cluster variance 
in galaxy properties at a fixed mass. Clearly larger sam- 
ples of well-studied clusters are needed at z > 0.5 over a 
range of cluster mass. 

5.3. The Integrated Luminosity of the Red Sequence 

In addition to determining the evolution (or lack 
thereof) in the shape of the luminosity function, a com- 



plete description of the growth of the red sequence 
also requires a measurement of how the total red se- 
quence light in clusters increases over time. To do this 
we integrate the observed red sequence LF for XMM- 
LSS J02182-05102 {Ltot,Rs) and compare it to that for 
clusters at 0.4 < z < 0.8 from EDisCS and at z = from 
SDSS. The results are shown in Figure [SJ Our LF depths 
at all redshifts have been constructed to have identical 
faint limits when corrected for passive luminosity evo- 
lution - indeed this drove our HAWK-I exposure times. 
For that reason, integrating the observed LFs yields an 
Ltot.RS measurement that extends to a roughly constant 
stellar mass limit. 

It is important to compare this cluster with its likely 
descendants at low reshift as the most massive clusters 
also have the red sequences with the highest Ltot.R.s (see 
Appendix). As described in fJHand as shown in Figure [3] 
the most hkely descendants of XMM-LSS J02182-05102 
are typical clusters at z < 1 with log(Af/MQ) '^ 14.3—15. 
In Figure [5] we have highlighted which of the EDisCS 
clusters are the likely descendants based on the projected 
mass evolution. 

We compare our clusters within a constant physical 
aperture of r = 0.75 Mpc. This does not take into ac- 
count that cl usters may preferent ially grow from the in- 
side out (e.g. iBalogh et al.ll2000l ). On the other hand, 
given the very uncertain mass of XMM-LSS J02182- 
05102 and the poor mass determination of some of our 
SDSS clusters - due to a small number of available spec- 
troscopic redshifts - using a constant metric aperture is 
more robust than one which scales as mass, e.g. i?200- 

In the left panel of Figure [S] it is clear that the Ltot,RS 
of the most luminous cluster red sequences at z < 1.62 
is comparable to within a few tenths of a dex. How- 
ever, we must account for the expected fading of the 
stellar populations in the red sequence galaxies as they 
age. To this end, we fade all of the cluster red sequences 
to z = using a simple stellar population formed at 
z = 2.35, which is consistent with the color evolution of 
the bri ght red sequence gala xies over this whole redshift 
range ("P apovich et al.l [20100 . The faded total luminosi- 
ties are shown in the right panel of Figure [5] Under the 
assumption that all the red sequence galaxies at all red- 
shifts have the same SFH, this roughly converts Ltot,RS 
into a stellar mass content on the red sequence-'^'^. We 
speculate on the evolution of XMM-LSS J02182-05102, 
by performing a least squares fit to the Ltot,RS values for 
XMM-LSS J02182-05102 and the two most 'massive (and 
most luminous) EDisCS clusters. This is shown as the 
dashed line in the right panel of Figure [S] 

Given the likely evolutionary path, the red sequence 
in XMM-LSS J02182-05102 grew by a factor of ~ 2 in 
light or stellar mass during the sa 4Gyr between z = 1.62 
and z = 0.6. In ^ we will discuss how to reconcile this 
rapid growth in the luminosity with the lack of shape 
evolution between z — 1.62 and z = 0.7. Assuming that 
the evolutionary path in Figure [5| continues to z = then 
the cluster will grow by an addition 50% between z = 0.6 
and z = or a factor of three in total since z = 1.62. 

As a cautionary note, given the rapid evolution in the 

^^ ISkelton et al.l | |2011I '| have pointed out that this assumption 
may be inappropriate given the heirarchical build-up of red se- 
quence galaxies but we assume it here for simplicity. 
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Figure 5. The total amount of light on the red sequence, Ltot,RSy in clusters at ^ < 1.62. The shaded histogram represents the 
distribution total red sequence luminosities from our local SDSS sample. In both panels, the triangles indicate the EDisCS clusters, with 
the larger triangles showing those clusters that are the likely intermediate redshift descendants of XMM-LSS J02182-05102 based on the 
likely mass growth from Figure|3] The errorbars on the EDisCS and XMM-LSS J02182-05102 point account for the statistical errors in the 
individual cluster LF determinations. Left panel: Lt^t^RSy computed by integrating the measured red sequence luminosity functions. Right 
panel: Both the EDisCS clusters and XMM-LSS J02182-05102 have had their Ltot,RS adjusted by fading the light on the red sequence to 
2 ~ by an amount expected for a simple stellar population with zj^^^ = 2.35. Using this M/L correction, we plot the effective stellar 
mass for galaxies at 2 = on the y-axis. The dashed line is a least squares fit to the Ltot,RS values for XMM-LSS J02182-05102 and its 
likely intermediate redshift descendants. 



shape of the LF at z < 0.6 it is likely incorrect to as- 
sume a constant amount of fading for all galaxies, as 
those added more recently to the red sequence will fade 
more rapidly vifith time. To approximate this effect we 
varied the amount of fading assuming a zjorm for XMM- 
LSS J02182-05102 ranging from z = 2 to 3 and for the 
EDisCS clusters ranging from z = 1 to 2. While the 
exact trend of Ltot,RS with z depends on the exact mix 
of Zform-, the change in maximum growth in Ltot,RS was 
~ 0.1 dex over the full redshift range. 

6. DISCUSSION 

Our first result, that there is a deficit of faint red galax- 
ies in XMM-LSS J02182-05102 may resuh naturally from 
a scenario in which galaxy star formation is quenched 
once a galaxy's total mass moves above IO^^Mq and it 
forms a hot gas halo. This naturally predicts that the 
most massive galaxies are quenched first and that there 
should b e relatively few low-mass passive galaxies at high 
redshift (|Gabor fc D ave 2012). 

In explaining our other findings on the growth of the 
red sequence over time we must explain some apparently 
contradictory results. First, we find that the LF of red 
sequence galaxies in clusters evolves very little in shape 
from z = 1.62 to z = 0.6. At redshifts lower than this, 
however, the shape evolves rapidly, such that the faint 
end slope becomes shallower, eventually matching the 
2 = value (Figure S]). At the same time. Figure [S] 
alone implies that the total stellar mass on the red se- 
quence appears to increase by a factor of ~ 2 during the 
~ 4 Gyr from z = 1.62 to z = 0.6 and then grows by 
only 50% over the remaining ~ 6 Gyr to the present day. 



In I Rudnick et al.l (j2009f) it was shown that adding the 
minimal possible number of galaxies to the red sequence 
at z < 0.8 would actually cause the predicted Ltot,RS 
in SDSS clusters to be too high. Resolving this discrep- 
ancy can be accomplished by assuming that a significant 
fraction of the stars in galaxies that are added to the 
red sequence at z < 0.8 end up as intracluster stars. At 
z > 0.7 we appear to have the opposite problem. The 
total light in the clusters must grow rapidly, but without 
changing the shape. 

At first glance, resolving this problem is difficult. 
There are no blue galaxies in the field at this redshift 
luminous enough to fad e onto the bright end of the red 
sequence in this cluster ()Papovich et al.|[2010[ ) , which im- 
plies that most of the mass must be added in the form 
of fainter galaxies. However, this would result in a shape 
change of the LF that is not seen between z = 1.62 and 
z = 0.7. A hint perhaps lies in the lack of red galaxies in 
XMM-LSS J02182-05102 brighter than Mi = -24.2, or 
Mj, fad e = —22.25, as has also been seen in iHilton et al.l 
(jlOOP). This is in contrast to the individual EDisCS 
clusters, which are populated with galaxies to at least 
0.5 magnitudes brighter than this, when accounting for 
the passive fadin g of the stellar popu lations (See indi- 
vidual LFs from IRudnick et al.l [20091 ). Therefore, de- 
spite the lack of blue luminous galaxies that could fade 
onto the red sequence, the bright red sequence popu- 
lation must grow between z — 1.62 and 0.7. One so- 
lution to this problem is for galaxies in the blue cloud 
at fainter magnitudes (~L*) to rapidly form stars and 
then be quenched, all on short enough timescales to pre- 
vent them from showing up in the observed population. 
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Figure 6. A demonstration of how merging and accretion will affect the evolution in the red sequence LF between z = 1.62 and 2 = 0.6. 
In all panels the vertical scaling of the luminosity functions are accurate in an absolute sense. The x-axis refers to the magnitude that 
galaxies are predicted to have at 2 = 0.6 assuming the SFHs given in the text. The red circles show the LF of XMM-LSS J02182-05102. 
The gray triangles show the composite EDisCS red sequence LF scaled in luminosity to the mean Ltat,RS of the likely descendant clusters 
at 2 0.6. These points are the same in each panel. The panels show a set of model predictions for the LF at 2 = 0.6 that assume that 
the cluster accretes enough recently quenched galaxies to account for the evolution in total red sequence light at Mi{z = 0.6) < —19. The 
squares in panel (a) show the predicted LF in the absence of merging of passive galaxies. The filled bands in panels (b), (c), and (d), 
represent how the predicted LF will change if every passive galaxy randomly merges two, three, or four times between 2 = 1.62 and 2 = 0.6. 
The extent of the bands corresponds to the 25 and 75% confidence intervals on the predicted LF resulting from a monte carlo simulation 
as described in the text. The extent of the bands in magnitude correspond to where the simulated LFs contain a median of 1 galaxy. The 
cluster LF shape and normalization at 2 = 0.6 can be approximately reproduced if every galaxy merges approximately three to four times 
over the intervening ~ 4 Gyr between 2 = 1.62 and 2 = 0.6. 



The amount of stars that need to be formed in such an 
episode, however, would require SFRs much in excess of 
what is measured for ga laxies in this field at z ~ 1.6 from 
MIPS 24/xm data (Tran et al.ll2010t ). 

6.1. The importance of mergers 

A more likely explanation is that the fainter blue galax- 
ies fall into the cluster, have their star formation sup- 
pressed somwhere during the process (in groups or the 
cluster itself), and migrate to the faint end of the red 
sequence, where they subsequently merge with other red 
sequence galaxies and increase their mass. It is reason- 



able that the SFRs of infalling galaxies are suppressed as 
they fall i nto the cluster envir onment, via ram-pressure 
stripping (Gunn & Gott 19 72|) or galaxy s t rangu lation 
(jLarson et al. 19801 . Indeed iPierre et all (|2011l ) have 
found that XMM-LSS J02182-05102 has a diffuse x-ray 
component indicating the presence of a heated (but pos- 
sibly underluminous) IGM. Such a scenario would have 
the effect of offsetting the addition of galaxies to the red 
sequence at faint magnitudes by merging them up the 
red sequence and hence preserving the shape of the LF, 
while at the same time allowing the cluster to grow its 
total red sequence stellar mass to agree with the most 
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massive EDisCS clusters. 

To test this we perform a simulation where we add 
galaxies to the red sequence between z = 1.62 and 
z = 0.6 and then test how the resultant red sequence LF 
shape and normalization depends on the different num- 
ber of mergers per galaxy. Because the shape and norm- 
lization are important for this test, we compare the LF 
for XMM-LSS J02182-05102 to the composite EDisCS 
LF normalized to the mean Ltot,RS of the likely descen- 
dant clusters. 

We assume that infalling field galaxies have a con- 
stant SFR until a time t^ut before they are rapidly 
quenched. This quenching happens a time tdeiay be- 
fore they are added to the cluster. This insures that 
they are red by the time they enter the cluster. During 
the quenching process they lose Ai magnitudes of inter- 
nal extinction. For the purpose of this test we choose 
tcut = 3 Gyr, tdeiay = 1-0 Gyr and Aj = . 7 maj 
(corresponding to Ay = 1 for a iCalzetti et al.l ()200^ 
attenuation curve), although we note that our results 
are not sensitive to the exact values adopted (see be- 
low). This scenario assumes that the SFRs of galax- 
ies are suppressed before falling into the cluster en- 
vironment, and that in the transformation to red se- 
quence galaxies they lose some amount of dust extinc- 
tion. This general scenario is consistent with previous 
works (e.g. JPog gianti ct al.' '2006': 'Dc Lucia et al.l 120071: 
Rudnick et al.l 12009: Dc Lucia ct al. 2007; McGee et ail 
2011^ and results in galaxies with colors on the red se- 
quence by z = 0.6. 

Over the 3.8 Gyr between z — 1.62 and 0.6 the clus- 
ter red sequence of XMM-LSS J02182-05102 has to in- 
crease its total stellar mass by a factor of 2 to match 
Ltot,RS for the likely descendants at z ~ 0.6. We split 
the necessary mass increase equally between 3 identical 
intervals of time over this redshift range. In each interval 
we d raw galaxies from the e volving blue field-galaxy LF 
from iSalimbeni et al.l (|2008D and evolve them using the 
SFH parameters above to predict their luminosities at 
z = 0.6. The Salimbeni LF is computed in the rest-frame 
B-band. We convert it to the rest-frame i-band assuming 
a constant SFH with 1 mag of extinction. The result of 
this process is shown in Figure [B^ as the yellow squares, 
which are compared to the XMM-LSS J02182-05102 and 
EDisCS LF. There we can see that such a scenario of 
simply accreting and quenching galaxies from the star- 
forming field population to make up the observed total 
luminosity increase would produce a faint end slope much 
steeper than what is observed in clusters at z = 0.6. 

We then merge each galaxy with a random second 
galaxy and perform this up to four times. We perform 
100 monte-carlo iterations of these mergers and show the 
25 and 75% limits of the distribution of the resultant LF 
as the yellow bands in Figure [6)3, c, and d for two, three, 
and four mergers per galaxy respectively. It appears that 
merging each red sequence galaxy three to four times 
with another red sequence galaxy results in an LF that 
is similar to that observed at z = 0.6. This implies a 
merger rate of ^ 1 per Gyr over this time period, inte- 
grated over all merger mass ratios. The faint-end slope 
for this model has the necessary turnover but the bright 
end is still underpopulated. If galaxies only merge twice 
then the faint-end slope is too steep. If galaxies merge 
four times, then the faint end turns over a little bit too 



rapidly but the bright end is somewhat above the obser- 
vations yet does not extend to the brightest magnitudes 
seen in EDisCS. 

These remaining challenges may imply that our model 
is too simplistic in its treatement of the mass ratios of 
mergers. For example, in a perfect world the mass ra- 
tios should be drawn from the LFs of infalling galax- 
ies and not from the general field population. Also, 
iDe Lucia et al.l (|2012f ) find that major mergers of mas- 
sive galaxies in clusters made be preferred based on their 
results using N-body simulations with a semi-analytic 
galaxy formation model. An additonal potentially im- 
portant unmodeled ingredient could the inclusion of the 
rapid buildup of stellar mass in massive cluster galaxies 
that are undergoing significant star formation but which 
must cease their star formation soon after we observe 
them (jTran et al.ll2010t ). Since the number of dusty star- 
forming galaxies increases rapidly towards high er red- 
shift, bot h in the field (e.g.UTe Floc'h et al.ll2005D and in 
clusters ("Saintonge et all 120081: iFinn et al.ll2010[ ). their 
rapid quenching may be a plausible mechanism to build 
up the bright end of the red sequence luminosity func- 
tion. 

It is worth noting that the exact number of mergers 
needed to match the LF is affected greatly by the slope 
of the faint-end of the blue field LF. If we assum e a = 
— 1.0 instead of -1.39 as in ISalimbeni et all ()2008D then 
the required number of mergers goes down to 1-2. The 
results of this model are also dependent on the faint limit 
from which we draw field galaxies. We use a faint limit of 
-14 but note that going to -16 would require two mergers 
to sufficiently deplete the faint-end of the LF. The results 
are less dependent on the SFH parameters tcut , tdeiay , 
and Ai. For example tcut < 2 Gyr and Ai < 0.5 both 
tend to favor four mergers, with three being ruled out.^^. 

The importance of merg ing in gr o wing the re d se- 
quenc e is consistent with iPapoviclil (|2011l ) and iLotd 
(|2011f ) who study this cluster with_WFC3 imaging from 
the CANDELS program. JLoti (pOll found that XMM- 
LSS J02182-05102 has a ~ 10 times higher merger rate 
than analogously selected galaxies in the coeval field, and 
that the close pairs are dominated by galaxies with a high 
mass ratio, i.e. minor mergers. The merger rate cited by 
iLotd (|2011f ) is 2 per Gyr, which is consistent with our es- 
timate of 1 per Gyr given the considerable uncertainties 
in both estimates, e.g. in the dependence of our derived 
mer ger rate on the field luminosity function. In addi- 
tion EaEQvicS l|2011f ) found using the same CANDELS 
data that the massive red sequence galaxies in XMM- 
LSS J02182-05102 require extensive dissipationless, and 
perhaps minor, mergers in order to simultaneously match 
the stellar mass, size, ellipticity, and color of cluster red 
sequence galaxies at z < 1. 

There is some direct evidence for merging of red se- 
quence galaxies in a merging cluster system at z — 0.8 
(Ivan Dokkum et al.l Il999t iTran et al.l 120051) . Likewise , 
iWhite et al.l (|2007l ) found that to evolve the Halo Oc- 
cupation Distribution for red luminous satellite galaxies 
(i.e. cluster red sequence galaxies) from z — 0.9 to 0.5 
requires that 1/3 of these galaxies must merge or un- 
dergo disruption in massive halos. This is again broadly 
consistent with the scenario that we propose here as the 

^* Ai = 0.7 corresponds to the typical opacity of galactic disks. 
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iWhite et alj (|2007D study includes redshifts where the 
growth in the red sequence has become more gradual. 
While merging is unlikely in a relaxed massive cluster 
due to the high velocities, XMM-LSS J02182-05102 will 
be rapidly growing with time, accreting other clusters 
and groups along the way. It is not relaxed, which im- 
plies that the galaxy-galaxy velocities may be lower than 
expected given the cluster mass. The relative veloci- 
ties in the infalling groups will also be lower than in 
the final cluster and at group scales the m erging cross- 
section may be quite high. Indeed iTran et al.l (|2008l ) 
find that the massive red galaxies in a set of merging 
groups at z ~ 0.37 are themselves experiencing dissipa- 
tionless mergers that will significantly grow their mass. 
By z = 0.8 the most massive clusters i n EDisCS hav e 
high velocity dispersions (Milvang-Jcnse ^et al.l 20081 ) . 
detected intracluster light (Gucnnou et al. 2011), and 
strong gravitational lensing (White ct al. 20051 . imply- 
ing that the cross section for merging may be much lower. 

6.2. Additional Caveats 

Our main uncertainty results from this study being 
based on one high redshift cluster. As pointed out in ij5.2l 
there are differences in the LF of massive clusters at 0.8 < 
z < 1.1 and we can expect that these differences might be 
more extreme in the less-developed cluster population at 
z > 1.5. Obviously, having larger samples of high redshift 
clusters with deep NIR data will be crucial if we wish 
to obtain a representative picture of galaxy evolution in 
clusters at these epochs. 

An additional uncertainty in our result regards the 
identification of red sequence galaxies with passively 
evolving ones. At z — 0.8 this is true at the ~ 85% 
level (Rudnick et al. in prep.). At z = 1.62. lTran et al.l 
(J2010) found that some of the red sequence galaxies were 
24/im emitters, indicating that a sma ll number of red 
gal axies could be dusty star formers. iPapovicbl (|2011[ ) 
and lOuadri et al.) (|2012D used combinations of rest-frame 
optical/NIR colors to separate red galaxies that are qui- 
esci ent from those that are obscured star-formers (see 
also I Williams et alll2009f )) and found that the contam- 
ination on the red sequence from dust-obscured objects 
is ^ 20% for bright red galaxies. Given the similar frac- 
tions of obscured red galaxies at both redshifts we con- 
clude that this will not cause a significant error in our 
conclusions. 

The mass estimate for our cluster is also uncertain, 
with a range between measurements of M ^ 7.7 x 10^"^ — 
1 X 10^'^M.Q. We adopted to lower of these two - from 
the x-ray detection - for the analysis in this paper. If 
we instead chose the mean of these two masses, then the 
projected mass growth tracks (yellow band in Figure [3]) 
would shift upward but remain roughly parallel to the 
original. The likely descendant clusters at z ~ 0.6 would 
then be the most massive clusters at every epoch with 
M > IO^^Mq although in reality such extreme descen- 
dants are unlikely due to the steepness of the cluster 
mass function. The EDisCS LF shape is not dependent 
on cluster mass (Rudnick ct_aL 2009), and so the ob- 
served lack of evolution in the shape should remain un- 
changed regardless of the expected descendants. If we 
adopt the mean mass then the expected luminosity evo- 
lution in i )5.3l and Figure [5] would be a factor 6 increase 
to z r^ 0.6 and almost no luminosity evolution to z = 0. 



As this is very extreme we are further confident that the 
lower mass estimate is more appropriate. 

A final caveat is our reliance on statistical background 
subtraction to determine membership. These should 
have equal or greater reliability when compared to pho- 
tometric redshifts ( §5.ip but spectroscopic redshifts are 
indeed sparse and will help to improve the constraints 
on the LF, especially at the bright end where continuum 
redshifts are feasible. Ample spectroscopy will also allow 
us to better measure the velocity dispersion and dynam- 
ical state of the cluster. Likewise, deep medium-band 
NIR imaging with 6 or 8-meter telescopes, can provide 
very precise photometric redshifts that would improve 
our analysis (jWhitaker et al.|[2011[ ). 

7. SUMMARY AND CONCLUSION 

We have presented the rest-frame i-band luminosity 
function of red sequence galaxies of XMM-LSS J02182- 
05102, a cluster at z = 1.62, as measured using deep 
HAWK-I observations from the VLT and a Ks se- 
lected catalog. At this redshift Y ~ J straddles the 
Balmer/4000A break and Kg is near in wavelength to 
rest-frame i. Our conclusions are as follows: 

• XMM-LSS J02 182-05 102 has a strong red sequence 
at bright magnitudes. Starting well brighter than 
our 90% completeness limit we find a lack of faint 
red galaxies. This conclusion is not dependent 
on our membership identification scheme as these 
galaxies simply do not exist in this area of the sky. 
The cluster red sequence hosts the objects with the 
reddest observed colors in this field and there are 
very few at faint magnitudes. 

• We derive a luminosity function for the red se- 
quence cluster members and compare it to anal- 
ogously constructed LFs for clusters at z = from 
SDSS and 0.4 < z < 0.8 from EDisCS. When 
corrected for passive evolution in the luminosities, 
we find that the shape of the XMM-LSS J02182- 
05102 LF is indistinguishable from the cluster red 
sequence LF at z = 0.7 and exhibits the same 
turnover to faint magnitudes that has been noted 
by other authors in z < 1 clusters. At z < 0.7 the 
faint end of the luminosity function starts to rise 
and has a near flat slope at z = 0. We also find 
that XMM-LSS J02182-05102 has a lack of lumi- 
nous red galaxies when compared to clusters from 
EDisCS. 

• The integral of the measured red sequence LF 
shows that LtotMS for XMM-LSS J02182-05102 is 
as high as the most luminous EDisCS cluster at 
0.4 < z < 0.8 and as the most luminous clusters 
in SDSS. However, when the expected evolution in 
the the stellar mass-to-light ratio is accounted for, 
XMM-LSS J02182-05102 has at least 3 times less 
stellar mass on the red sequence compared to its 
likely descendant clusters in SDSS and a factor of 
2 less stellar mass than its likely descendants in 
EDisCS. 

• We attempt to explain the large growth in the to- 
tal luminosity between z = 1.62 and z = 0.6 while 
simultaneously preventing the buildup of faint red 
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galaxies and also growing massive red ones. We test 
a simple model in which the cluster accretes galax- 
ies from the blue field population that are quenched 
and then merge on the red sequence. If every clus- 
ter red sequence galaxy mergers ^3 — 4 times 
with another red sequence galaxy, we can satisfy 
all of the above constraints. This corresponds to a 
merger rate of ~ 1 Gyr~^ integrated over all mass 
ratios. This also agrees with independent merger 
rates from this cluster. 

There are some limitations to our analysis. First, our 
membership information is limited due to the difficulty 
in obtaining spectroscopic membership information. NIR 
spectrographs on 6 and 8-meter telescopes will be crucial 
here. We also only have one cluster with an LF of mod- 
erate signal-to-noise. Obviously, increasing the sample 
of clusters will allow for a more precise determination of 
the mean cluster LF at these redshifts and will let us de- 
termine if XMM-LSS J02182-05102 is typical of the high 
redshift cluster population. 
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Figure Al. A comparison of the dynamical and weak-lensing mass with the integral of the red sequence luminosity L^ot RS- The circle 
and triangle points are values for the EDisCS clusters 0.4 < z < 0.8. Left panel: The EDisCS masses are derived from weak lensing 
estimates (Clowe ct al. 2006), where the blue points are those where substructures may contaminate the lensing signal. The horizontal 
dashed line and yellow band (left panel) and red point (right panel) gives Ltot,RS and its uncertainty for XMM-LSS J02182-05102, where 
the luminosities have been passively faded to z = 0.6. Right pa nel: The EDisC S masses in this panel have been derived from the velocity 
dispersion using 30-50 memb ers per system fro m IMil vang- Jen sen et aLl l|2008l). The red solid square represents the x-ray derived mass 
for XMM-LSS J02182-05102 (Pi erre et al.ll20lID . The open square connected to the solid by a dotted line indicates the dynamical mass 
estimate from Papovich ct al. ( 2 0101 ). The disagreement between the two likely reflects the unrelaxed dynamical state of the system. All 
luminosities are computed in a circular aperture with a radius of 0.75 physical Mpc. The EDisCS luminosities have all had a small passive 
evolution correction to z = 0.6. The correlation of Ltot,RS a-^d cluster mass is significant for lensing mass, but less so for mass from the 
velocity dispersions. 

APPENDIX 
THE USE OF THE INTEGRATED RED SEQUENCE LIGHT AS AN ESTIMATOR OF THE DYNAMICAL MASS. 

Since t he dy namical state of our cluster is likely not relaxed, we seek to estimate the mass using a different method . 
In Figure lATJ we show the correlation of Ltot.R s in a r = 0.75 Mpc apertur e vs. weak lensing mass (jClowe et al.ll2006l ) 
and mass derived from the velocity disperion (jMilvang- Jensen et al.| [2008[) for the EDisCS clusters at 0.4 < z < 0.8. 
There is a clear correlation between lensing mass and Ltot,RS but a significantly poorer relation with mass derived 
from the velocity disperion, even though this dispersion is measured from 30-50 members. The horizontal line shows 
the value for XMM-LSS J02182-05102, where we have corrected the luminosities for passive evolution to z = 0.6. If 
we assume the same ratio of dark matter mass to stellar mass on the red sequence at z = 0.6 and z — 1.62 then it 
would appear that XMM-LS S J02182-05102 has M ~ 2. 5 x 10^"^ Mf^ which is in excellent agreement with the value 
of M ~ 1 — 4 X 10^** M0 from iPapovich et al.l (|2010f) and lPierre et al.l (|2011f) as determined from the x-ray luminosity 
on the low side and the very uncertain velocity dispersion (ct = 860 ± 490km s"'^) on this high side. This opens the 
possibily of using Ltot,RS as a proxy for mass in clusters that we are now selecting with this method, but which do 
not have spectroscopy sufficient to determine a velocity dispersion. 



